147 research outputs found
A survey of partial differential equations in geometric design
YesComputer aided geometric design is an area
where the improvement of surface generation techniques
is an everlasting demand since faster and more accurate
geometric models are required. Traditional methods
for generating surfaces were initially mainly based
upon interpolation algorithms. Recently, partial differential
equations (PDE) were introduced as a valuable
tool for geometric modelling since they offer a number
of features from which these areas can benefit. This work
summarises the uses given to PDE surfaces as a surface
generation technique togethe
A Moving Boundary Flux Stabilization Method for Cartesian Cut-Cell Grids using Directional Operator Splitting
An explicit moving boundary method for the numerical solution of
time-dependent hyperbolic conservation laws on grids produced by the
intersection of complex geometries with a regular Cartesian grid is presented.
As it employs directional operator splitting, implementation of the scheme is
rather straightforward. Extending the method for static walls from Klein et
al., Phil. Trans. Roy. Soc., A367, no. 1907, 4559-4575 (2009), the scheme
calculates fluxes needed for a conservative update of the near-wall cut-cells
as linear combinations of standard fluxes from a one-dimensional extended
stencil. Here the standard fluxes are those obtained without regard to the
small sub-cell problem, and the linear combination weights involve detailed
information regarding the cut-cell geometry. This linear combination of
standard fluxes stabilizes the updates such that the time-step yielding
marginal stability for arbitrarily small cut-cells is of the same order as that
for regular cells. Moreover, it renders the approach compatible with a wide
range of existing numerical flux-approximation methods. The scheme is extended
here to time dependent rigid boundaries by reformulating the linear combination
weights of the stabilizing flux stencil to account for the time dependence of
cut-cell volume and interface area fractions. The two-dimensional tests
discussed include advection in a channel oriented at an oblique angle to the
Cartesian computational mesh, cylinders with circular and triangular
cross-section passing through a stationary shock wave, a piston moving through
an open-ended shock tube, and the flow around an oscillating NACA 0012 aerofoil
profile.Comment: 30 pages, 27 figures, 3 table
Anomalous wave structure in magnetized materials described by non-convex equations of state
Agraïments: Institute for Pure and Applied Mathematics (UCLA) 2012 program on "Computational Methods in High Energy Density Plasmas.We analyze the anomalous wave structure appearing in flow dynamics under the influence of magnetic field in materials described by non-ideal equations of state. We consider the system of magnetohydrodynamics equations closed by a general equation of state (EOS) and propose a complete spectral decomposition of the fluxes that allows us to derive an expression of the nonlinearity factor as the mathematical tool to determine the nature of the wave phenomena. We prove that the possible formation of non-classical wave structure is determined by both the thermodynamic properties of the material and the magnetic field as well as its possible rotation. We demonstrate that phase transitions induced by material properties do not necessarily imply the loss of genuine nonlinearity of the wavefields as is the case in classical hydrodynamics. The analytical expression of the nonlinearity factor allows us to determine the specific amount of magnetic field necessary to prevent formation of complex structure induced by phase transition in the material. We illustrate our analytical approach by considering two non-convex EOS that exhibit phase transitions and anomalous behavior in the evolution. We present numerical experiments validating the analysis performed through a set of one-dimensional Riemann problems. In the examples we show how to determine the appropriate amount of magnetic field in the initial conditions of the Riemann problem to transform a thermodynamic composite wave into a simple nonlinear wave
A new ghost cell/level set method for moving boundary problems:application to tumor growth
In this paper, we present a ghost cell/level set method for the evolution of interfaces whose normal velocity depend upon the solutions of linear and nonlinear quasi-steady reaction-diffusion equations with curvature-dependent boundary conditions. Our technique includes a ghost cell method that accurately discretizes normal derivative jump boundary conditions without smearing jumps in the tangential derivative; a new iterative method for solving linear and nonlinear quasi-steady reaction-diffusion equations; an adaptive discretization to compute the curvature and normal vectors; and a new discrete approximation to the Heaviside function. We present numerical examples that demonstrate better than 1.5-order convergence for problems where traditional ghost cell methods either fail to converge or attain at best sub-linear accuracy. We apply our techniques to a model of tumor growth in complex, heterogeneous tissues that consists of a nonlinear nutrient equation and a pressure equation with geometry-dependent jump boundary conditions. We simulate the growth of glioblastoma (an aggressive brain tumor) into a large, 1 cm square of brain tissue that includes heterogeneous nutrient delivery and varied biomechanical characteristics (white matter, gray matter, cerebrospinal fluid, and bone), and we observe growth morphologies that are highly dependent upon the variations of the tissue characteristics—an effect observed in real tumor growth
- …